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ABSTRACT 

The subject of this paper is the derivation of the integrated Sachs-Wolfe (iSW) effect in cos- 
mologies with coupled dark matter and dark energy fluids. These couplings influence the iSW- 
effect in three ways: The Hubble function assumes a different scaling, the structure growth 
rate shows a different time evolution, and in addition, the Poisson equation, which relates 
the density perturbations to fluctuations in the gravitational potential, is changed, due to the 
violation of the scaling p oc a~ 3 of the matter density p with scale factor a. Expemplarily, I 
derive the iSW-spectra for a model in which dark matter decays into dark energy, investigate 
the influence of the dark matter decay rate and the dark energy equation of state on the iSW- 
signal, and discuss the analogies for gravitational lensing. Quite generally iSW-measurements 
should reach similar accuracy in determining the dark energy equation of state parameter and 
the coupling constant. 
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1 INTRODUCTION 



Recently, cosmological models with coupling terms in the evo- 
lution equations for the dark matter and the dark energy den- 
sity have attracted interest, because these couplings alleviate the 
coincidence problem in an elegant way dChimento etall 120031; 



^ ■ l Olivares et alj2007l;|Pettorino & Baccigalupill2008l;lBoehmer et al.l 



2008; He & Wandl2008l) . 

The integrated Sachs - Wolfe (iSW) effect dSachs & Wolfd 
19671: iRees & Sciamal 19681: iHu & Sugivamalll994 : ICooravll2002l ; 
Schafer & Bartelmannl 120061) , which refers to the frequency 



change of cosmic microwave background (CMB) photons if 
they cross time evolving gravitational potentials, is a direct 
probe of da rk energy because i t van ishes in cosmologies with 
£2 m = 1 dCrittenden & Turokl 1 19961) . By now, it has been 
detected with high significance with a number of different 
tracer objects |Boughn & C rittenden 120031: IFosalba et alj 12003 



Giann antonio et all 
20081) . 



2006; Rassatetal. 2007; Giannantonio et al 



In this paper, I would like to focus on the derivation of the 
integrated Sachs-Wolfe effect for cosmological models with cou- 
pled dark matter and dark energy fluids and to show that iSW-effect 
acquires a dependence on the Hubble function, the matter density 
parameter, the growth rate and their derivatives. In this respect, the 
iSW-effect in coupled models is much richer than in models with 
non-interacting constituents, where the iSW-effect simply measures 
the line of sight integrated derivative d(D + /a)/da of the growth 
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function D + (a). Physically, the changes in the iSW-formulae are 
due to the fact that the matter density p„, does not develop oc a~ 3 in 
coupled models, and the relation £2 m (a)/£2 m = Hi/fPiaj/a 3 does 
not hold. A consequence of this is a different (time-evolving) rela- 
tion between the gravitational potential and the overdensity field; 
and it is not clea r how mechanism wa s incorporated in the deriva- 
tion presented bv lOlivares et al.l |2008). 

After extending the iSW-formulae to coupled cosmological 
models in Sect. [2] I compute the power spectrum of the iSW-effect 
in Sect.f3]for a phenomenological model in which dark matter de- 
cays into dark energy. A summary of my results is compiled in 
Sect. [4] I consider spatially flat homogeneous dark energy cos- 
mologies with adiabatic initial conditions in the cold dark matter 
field. Specific parameter choices are Hq = 100/? km/.s/Mpc with 
h = 0.72, £l m = 0.25, cr 8 = 0.8 and n s = 1. 



2 COSMOLOGY WITH COUPLED DARK FLUIDS 
2.1 Decay of dark matter into dark energy 

As an examp le, I compute th e iSW- spectra for the coupled model 
proposed by iBoehmer et alj d2008l) . in which cold dark matter 
(CDM) decays into dark energy: The evolution of the dark matter 
density p,„ and the dark energy density (with an equation of state 
parameter w) is described by the system of differential equations, 
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Table 1. Summary of the four dark energy models considered in this paper, 
in terms of the dark energy equation of state parameters wq, w a and the 
CDM decay constant T. 

with the cosmic time t as the time variable. The Friedmann con- 
straint is given by 3H 2 = p m + p$. The coupling constant T corre- 
sponds to the CDM decay rate in units of the Hubble-constant Hq. 
Models with stable CDM and hence uncoupled fluids are recovered 
by setting T = 0. 1 reformulate the time derivatives in the system of 
equations as derivatives with respect to the scale factor a for intro- 
ducing t he common parameter i sation for the dark ener gy equation 
of state dTurner & White! 19971 ; Kinder & Jenkinsll2003l) . 

w(a) = w + (1 - a)w a . (2) 

With solutions forp,„(a) and p^(a), the Hubble function H(a) is can 
be obtained by integrating, 

dH 2 f s 

-t— = ~ \Pm(a) + P4,{a) [1 + w(a)] J , (3) 

which in addition gives the definition of the critical density 
Pciit( a ) = 3H 2 (a)/(8nG), with Newton's constant G. Initial con- 
ditions are defined at the present epoch, H(a = 1) = Hq, £l m (a = 
1) = f2,„ and f2^(a = 1) = l-Q. m . The comoving distance^- follows 
from the solution for H(a), 

C 1 da 

X (d) = c -j— (4) 
J a a-H(a) 

c denoting the speed of light, and the density parameters can 
be defined by normalising Sl m (a) = p m (a)l p crlt (a) and fi^(a) = 
p,p{a)l p a \ t (a). TableQ]gives an overview over the dark energy mod- 
els considered in this paper, in terms of their dark energy properties 
w , w a and their CDM decay rate T. In the following I focus on 
spatially flat models with Cl„, + Q$ = 1, which is conserved in the 
time evolution of eqn. {TJ. 

2.2 Structure growth in dark energy cosmologies 

The homogeneous growth of the overdensity field, 6(x, a) = 
D + (a)S(x, 1) is described by the growt h function D+(a), which is 
a solution to the diffe rential equation dWang & Steinhardj |l998l ; 
iLinder & Jenkinsll2003h . 

^Ld + + -(3 + ^-^] ^-D + = ^-n m (a)D + (a). (5) 
aa a \ a ma j aa 2a^ 

In the standard cold dark matter (SCDM) cosmology with Q.„, = 1 
and 3 + d In H/d In a = | , the solution for D + (a) is simply the scale 
factor, D + (a) = a. Eqn. (J5} can be applied for describing structure 
growth in coupled models as well, because the overdensity field, 
being the ratio between the density perturbation and the mean den- 
sity, is not affected by e.g. CDM decay, as long as the dark energy 
fluid can be considered to be homogeneous and not to be influ- 
encing local structure formation. Nevertheless, coupled models in- 
fluence the growth equation by the scaling of the Hubble function 
H(a) and the time evolution of the density parameter Ci m (a). 




scale factor a 

Figure 1. The growth function D + (a) as a function of scale factor a, for 
ACDM (solid line), ArCDM with T = \ (dashed line), 0CDM (dash-dotted 
line), and ftfCDM with T = | (dotted line). In addition, the growth D + (a) = 
a for the SCDM cosmology is plotted (solid straight line). 

Additionally, the initial conditions of the growth equation are 
different, as the cosmology does not have an asymptotic SCDM 
phase with H oc a~ 3/2 at early times. The inital conditions can be 
specified by assuming an asymptotic power-law growth D + (a) = 
a", a > 0. Substitution into eqn. {5} yields a quadratic equation for 
a, whose positive solution specifies the initial conditi ons £> + (fl,) = 
qf an d dD + (o,)/da = ara" -1 at some early time a, (c.f. lSchafer et alj 
120081) . 

Fig. Q] shows growth functions for the four exemplary cos- 
mologies. Evolving dark energy has the property of suppressing 
structure formation at an earlier time, which can be partially com- 
pensated by CDM decay (because the gravitational fields gener- 
ated by the overdensity 6(x) are stronger if Sl„,(a) has a higher 
value), causing degeneracies between the equation of state parame- 
ters and the CDM decay rate: The ACDM and r CDM-models, for 
instance, are almost indistinguishable. 

It is worth noting that models with decaying CDM are gen- 
uinely different from dark energy models concerning structure 
growth and probes which use gravitational interaction such as grav- 
itational lensing or the iSW-effect. While it is always possible to 
construct an equation of state w(a) for a dark energy model with 
stable CDM that gives the identical Hubble function as a model 
with decaying CDM, this degeneracy is broken in the evolution of 
the density parameter Sl m (a) and the growth function D + (a) which 
both would be different in the two cases. Theoretically, combining 
observations of cosmic structure growth and the expansion history 
should make it to distinguish between the two families of mod- 
els, which geometrical probes such as supernova observations alone 
would not be able to do. 



2.3 iSW-effect in coupled models 

The iSW-effect is caused by gravitational interaction of a CMB 
photon with a time-evolving potential O. The frac tional perturba- 
tion t of the CMB temperature 7cmb is given bv dSachs&Wolfel 
fl967h 

t = ^ = -T \ dl — = — - dxa l H(a) — , (6) 

Tcmb c 1 J dr] c 3 Jo da 
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where 77 denotes the conformal time. In the last step, I have re- 
placed the integration variable by the comoving distance x, which 
is related to the conformal time by dx = —cdr) = —cdt/a, and the 
time derivative of the growth function has been rewritten in terms 
of the scale factor a, using the definition of the Hubble function 
da/At = aH(a), with the cosmic time t. The gravitational potential 
O follows from the Poisson equation in the comoving frame, 



AO = 4nGa-p,„(a) d 



(7) 



AO = -a 2 H 2 (a)Q„,(a) 6. 



(8) 



where Newton's constant G is replaced with the critical density 
Pcrit(fl) = 3H 2 (a)/(SnG), p m {d) = n m (a)p crit (a), 
3 

2 C 

Substitution yields a line of sight expression for the linear iSW- 
effect t (using the Born-approximation and assuming linear struc- 
ture formation), 

3//;; f ! '< 



Hi r XH . dQ . 

dx a H(a) — A S, 
r - Jo da 



(9) 



with the inverse Laplace operator A 1 solving for the potential. The 
function Q(a) is given by 



Q(a) = a'h (a)£l m (a)D + (a), 

with h(a) = H(a)/Ho, and has the derivative dQ/da, 



da all din a 



d in Q.„, d In D + 



din a 



din a 



(10) 



(11) 



i.e. the iSW-effect measures the sum of the scalings of the Hubble 
function, of the matter density and of the growth function, inte- 
grated along a line of sight. 

The time evolution of the iSW-source term Q(d) and its deriva- 
tive dQ/da are depicted in Fig.[2] All models intersect at Q(a) = D„, 
at the current epoch, and models with evolving dark energy or CDM 
decay show higher values for Q(a) in the past. Both the models with 
decay and the models with evolving dark energy assume similar 
values for dQ/da at low and high redshifts, respectively, and evolv- 
ing dark energy shifts the curves horizontally. Again, the model 
with evolving dark energy and CDM decay produces the highest 
amplitudes for the derivative dQ/da. 

Physically, the expression for Q(a) corresponds to the three 
ways in which models with CDM decay generate higher amplitudes 
for the iSW-effect, compared to models with stable CDM: Apart 
from the growth factor D+{a) and the Hubble function H(a) (due to 
the time derivative of D + ), the matter density has had higher values 
in the past, and hence the fluctuations in the gravitational potentials 
generated by the overdensity field S are stronger, which is described 
by Q.„,(a). Further interesting points include: 

• The derivative dQ/da in fact vanishes in the SCDM cosmol- 
ogy with Q,„(a) = 1, D + (a) = a and h(a) = a~ 3/2 , because Q(a) = 1 
in this case. It is perhaps even more illustrative to look directly 
at the expression for dQ/da (eqn. lilt which is zero in SCDM 
cosmologies because (i) the logarithmic derivative of the Hubble 
function H(a) is equal -3/2 (reflecting the matter domination), (ii) 
the logarithmic derivative of the matter density parameter vanishes 
(due to the stability of CDM), and the logarithmic derivative of 
the growth function D + (a) is equal to one (a combination of the 
two). From this point of view, both the growth proportional to a 
and stable CDM are the necessary conditions that make the linear 
iSW-effect disappear in Einstein-de Sitter universes. 

• The two cosmologies ACDM and r CDM, which are almost 
degenerate in terms of their growth rates, give rise to very different 
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Figure 2. Time evolution Q(a) of the iSW source term (thick lines) and its 
derivative dQ(a)/da (thin lines) as functions of scale factor a, for ACDM 
(solid line), A r CDM with r = I (dashed line), 0CDM (dash-dotted line), 



and 0rCDM with r = £ (dotted line). 



Q(a) and dQ/da functions, which is a result of the their dependence 
on H(a) and fi,„(a), and hence to a different iSW-effect. 

• In models with no coupling between the dark fluids, the matter 
density decreases p m (a) cc a~ 3 , and 



flm(a) 

Qui 



a 3 H 2 (a) 



(12) 



applies. This allows the reformulation of the comoving Poisson 
equation in the familiar form, 



3H 2 Q„, 
AO = — S, 

la 



(13) 



Substitution into the line of sight integral yields the familiar equa- 
tion for the iSW-effect in spatially flat dark energy cosmologies: 



c 3 Jo 



d X a 2 H(a) 



d D + 
da a 



A" 1 S, 



(14) 



and Q(a) reduces to £l ln D + (a)/a with a constant Q. m . 

• An analogous consideration applies to gravitational lensing: 
The weak lensing convergence k in coupled models reads 



3 r 



dx a 2 H 2 (a)Q. m (a) G(x)x D+S, 



(15) 



which replaces the common expression jSchneideretai]|l992l: 
iBartelmann & Sc hneider 200lb, 



Q t/2Q 

m 

2c 2 



d X G(x)x—6- 

a 



(16) 



G(x) is the lensing efficiency function of the background galaxy 
sample, which itself is influenced by the Hubble function H{a). 

The derivation above holds for general coupled cosmologies and is 
not restricted to the case of CDM decay, for which I exemplarily 
compute the iSW-spectra in the next chapter. 
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3 ISW ANGULAR POWER SPECTRA 

In summary, the line of sight integrals for the iSW-temperature per- 
turbation t and the galaxy density y read: 



r 

Jo 



d X a 2 H(a)-^ V , 
da 



o 

b d XP {z)— D + (x)6, 
Jo 



d X 



(17) 



(18) 



where I have defined the dimensionless potential if = h~ l 6ld 2 H , 
rescaled with the square of the Hubble distance d H = c/Hq for 
convenience. The weighting functions 



W T (x) 
W 7 (x) 



3 , dQ 
-a 2 H(a)-=-, 
c da 

p(z)^bD + (a), 
dx 



(19) 



(20) 



can be identified, which allow the expressions for the iSW-auto 
spectrum C TT (£), the iSW-cross spectrum C Ty {£) and the galaxy 
spectrum C ry( f) to be written in a compact notation, applying a 
Limber-projection dLimbe J 195 4*) in the flat-sky approximation, for 
simplicity: 



Crrd) = f 
JO 

- I" 

Jo 

- f 

Jo 



CryiO 



dx 



dx 



dx 



W 2 (x) 



X 

w T ( X W y { x ) 



P w (k = tlx), 



P<rs(k = llx\ 



x 2 



(21) 



(22) 



(23) 



with the cross-spectrum P^(k) = P S6 (k)/(d H k) and the spectrum 
Pipipik) = Pss(k)/(d H k) 4 of the potential if. 

For the CDM power spectru m I make the ansa tz P(k) oc 
k"'T 2 (k), with the transfer function teardeen et"aDll986h . 

ln(l + 2.34o) / 

where the wave vector q is given in units of the shape parameter 
Q. m h. P(k) is normalised to the value erg on the scale R = 8 Mpc/h, 



' + (I6.I4) 2 + (5A6qf + (6.71?) 4 ) ? 



2n 2 J 



dk k 2 W 1 (kR)P(k), 



(24) 



with a Fourier-transformed spherical top-hat W(x) = 3ji(x)/x as 
the filter function. j e (x) denotes the spherical Bes sel function of the 
first kind of order I dAbramowitz & Stegunll972h . The contribution 
of nonlinear structure formation to th e power spectrum P(k) was 
described with the model proposed by Smith et al. (2003), yielding 
the power spectrum Pss(k), where the parameterisation of the non- 
linear contribution to the fluctuation amplitude with Q.„,(a) as the 
time variable is suited for coupled models. 

I use the redshift distrib ution of the main galaxy sample of the 
Dark UNiverse ExploreiQ i Refregie r & the DU NE collab oration! 
2008), which will observe half of the sky out to redshifts of 
unity and which wi ll be of particular use for iSW-observations 
dPouspis et alj2008h : 



p(z)dz = p 



— I exp 



— I dz with 



Po P \P 



(25) 



with p = 3/2 and zo = 0.64, which results in a median redshift of 



1 http://www.dune-mission.net/ 
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Figure 3. iSW auto power spectra C TT (() as a function of inverse angular 
scale I, for ACDM (solid line), A r CDM with r = \ (dashed line), <pCDM 
(dash-dotted line), and r CDM with T = 5 (dotted line). 



Zmed = 0.9. For simplicity, the bias b is assumed to be non-evolving 
and equal to unity. 

The angular spectra C TT (£) and C Ty (€) are shown in Figs. [3] 
and [4] respectively, for the decaying CDM cosmology outlined in 
Sect. [2] Differences between the cosmologies considered are not 
strongly scale dependent as expected for changes on the homoge- 
neous level and observing linear structure formation, and amount to 
half an order of magnitude in the auto spectrum C TT (C). The cross 
spectrum C Ty {£) exhibits differences up to a factor of two. Both 
spectra seem to be equally sensitive to the dark energy equation of 
state as well as to the coupling term T, and these two parameters 
are naturally degenerate with the fluctuation amplitude erg of the 
density field. For the chosen redshift distribution, most of the cross 
signal originates at a redshift of z — 1.0 (a = 0.5, c.f. Fig.[2]l, where 
the source function dQ/da is significantly stronger in models with 
either evolving dark energy or CDM decay, compared to ACDM. 
The spectra show degeneracies between w, T and cr 8 , such that a 
precision measurement will have to rely on a good prior on erg in 
order to unlock the sensitivity of the iSW-effect on w and T. 



4 SUMMARY 

In this paper, I extend the expressions for the iSW-effect to cosmo- 
logical models with couplings between the dark matter and dark 
energy densities, and compute the auto and cross correlation angu- 
lar spectra for a coupled cosmology, in which dark matter decays 
into dark energy. 

• The iSW-effect reflects couplings between the dark matter 
and dark energy fluids by three mechanisms: The Hubble func- 
tion H(a) shows a different scaling in the past, the growth equa- 
tion has different initial conditions and a different time evolution, 
and the Poisson equation, which relates the gravitational potential 
to the overdensity field, is affected by the evolution of the mat- 
ter density Q.„,(a). The iSW-formula reduces to the familiar form if 
Q. m (a)/£l„, = HllH 2 (a)la i is substituted, which holds in uncoupled 
models. 

• The iSW-spectrum C TT (£) shows higher amplitudes in mod- 
els with couplings as well as in models with evolving dark energy, 
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Figure 4. iSW cross power spectra C Ty (t) as a function of inverse angular 
scale €, for ACDM (solid line), A r CDM with T = \ (dashed line), 0CDM 
(dash-dotted line), and r CDM with T = 5 (dotted line). 

compared to models with stable matter or those with a cosmolog- 
ical constant A. For the (extreme) choice of T and w the cross- 
spectrum C Ty (t) assumes amplitudes which are twice as large in 
dark energy models with evolving w or in coupled models, i.e. the 
sensitivity of the iSW-effect towards the mean dark energy equation 
of state parameter and coupling terms is comparable. The shape of 
the spectra suggests degeneracies between w, T and cr 8 , which all 
increase the overall power. For that reason, good independent priors 
on erg are essential to constraints on w and T from the iSW-effect. 

• Because of the weakness of the iSW-effect it will be difficult 
derive strong constraints on size of the interaction term, or the spe- 
cific type of the coupling. Gravitational lensing, on the contrary, is 
equally affected by a violation of the scaling of the matter density, 
and is likely to yield constraints comparable to those on the dark 
energy equation of state parameters, possibly to an accuracy of a 
few percent with upcoming surveys dLa Vacca & Colombo! 2008). 
A companion paper considers structure formation in coupled mod- 
els and examins the measurement of the CDM decay rat e T with 
weak lensing bispectrum tomography JSchafer et ai]|2008h . 
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